Odd-frequency Cooper pairs and zero-energy surface bound states in superfluid 3 He 
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We study the odd-frequency Cooper pairs formed near the surface of superfluid 3 He. The odd- 
frequency pair amplitude is closely related to the local density of states in the low energy limit. We 
derive a formula relating explicitly the two quantities. This formula holds for arbitrary boundary 
condition at the surface. We also present some numerical results on the surface odd- frequency pair 
amplitude in superfluid 3 He-B. Those analytical and numerical results allow one to interpret the 
midgap surface density of states, observed recently by transverse acoustic impedance measurements 
on superfluid 3 He-B, as the manifestation of the surface odd-frequency state. 
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I. INTRODUCTION 

The spin and orbital degrees of freedom of the Cooper 
pairs are customarily classified into two categories: spin- 
singlet even-parity and spin-triplet odd-parity. This clas- 
sification is based on the requirement that the equal-time 
pair amplitude is antisymmetric under simultaneous ex- 
change of spin and position variables. Berezinskii has 
proposed a concept of even- and odd-frequency symme- 
tries associated with the relative-time dependence of the 
pair amplitude^ The even- frequency state is character- 
ized by the pair amplitude symmetric in frequency, so 
that it can be subdivided into the same symmetry classes 
as in the equal-time case. On the other hand, the odd- 
frequency state can have the symmetries of spin-singlet 
odd-parity and spin-triplet even-parity, in contrast to 
the even- frequency state. The even-frequency symme- 
try class includes almost all the paring states so far dis- 
cussed. For instance, the spin-singlet s-wave state dis- 
cussed by Bardeen, Cooper, and Schricffer for supercon- 
ducting metals^ the spin-singlet d-wave state for cuprate 
superconductors (2, and the spin-triplet p-wave states for 
superfluid 3 He^ and for S1-2RUO4— belong to this cat- 
egory. The possibility of the odd-frequency pairing has 
been proposed in the context of cuprate superconductors^ 
and heavy Fermion superconductors^ At present, how- 
ever, the odd-frequency superconductor or superfluid has 
not yet been established experimentally. 

Bergeret, Volkov, and Efctov pointed out that the odd- 
frequency Cooper pairs are generated by spatial inho- 
mogeneity even without interaction responsible for odd- 
frequency pairing^ After that, several theoretical studies 
have been reported on the odd-frequency states in super- 
conducting j unctions An intriguing example of the 
junctions is a diffusive normal metal/spin-triplet super- 
conductor system. In this system, the proximity-induced 
superconductivity in the diffusive normal metal is dom- 
inated by the odd-frequency spin-triplet s-wave Cooper 
pairs because of isotropization by impurity scattering. 9 



Such an anomalous proximity effect is also expected in 
superfluid 3 He partly filled with aerogel*^ In general, 
there can coexist even- and odd-frequency Cooper pairs 
near interfaces or surfaces owing to broken translational 
symmetry. The inhomogeneous superconducting and su- 
perfluid systems therefore provide a good starting point 
to develop our understanding of the odd-frequency state. 

In this paper, we discuss the odd-frequency Cooper 
pairs generated near the surface of superfluid 3 He. The 
spin-triplet p-wave state of bulk superfluid 3 He is ar- 
guably the best understood non-s-wave stated It is well 
established that the non-s-wave states have pronounced 
sensitivity to surface scatteringJ^r— Unlike the s-wave 
state, the quasiparticle scattering at the surface causes 
substantial pair-breaking. As a result, the gap func- 
tion varies in space over several coherence lengths from 
the surface. In addition, the Andreev bound states are 
formed at the surface, yielding a characteristic low energy 
structure in the surface density of states below the bulk 
gap. The surface Andreev bound states have recently 
attracted renewed interest from the aspect of their Ma- 
jorana naturei^Ir— The purpose of the present work is to 
demonstrate the existence of the surface odd-frequency 
Cooper pairs in superfluid 3 He. Recently, Aoki et alr^ 
and Murakawa et a/. 20 ' 21 performed transverse acoustic 
impedance measurements to probe the surface of super- 
fluid 3 He and found an evidence of the low energy den- 
sity of states due to the formation of the surface An- 
dreev bound states in the B phase of superfluid 3 He. We 
show that the observed midgap density of states can also 
be interpreted as the manifestation of the surface odd- 
frequency Cooper pairs in superfluid 3 He-B. 



II. QUASICLASSICAL THEORY 

Our analysis is based on the quasiclassi- 
cal Green's function theory of superfluidity and 
superconductivity^ - — We apply it to a semi-infinite 
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system occupying the space z > 0. The surface at z — 
may have atomic-scale irregularities, though we assume 
it to be macroscopically flat. The surface roughness of 
atomic scale, which gives rise to diffuse quasiparticle 
scattering, is inevitable in actual systems. 

The quasiclassical Green's functions, g and /, are a 
function of (p, e, z), where p is a unit vector in the direc- 
tion of the Fermi momentum and e is a complex energy 
variable. One can express g and / in the form 28 
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where J-(p, e, z) and J-{p, e, z) are 2x2 spin-space ma- 
trices and the notation X denotes the transformation 
X(p, e,z) = X(— p, — e*, z)* . The spatial dependence of 
T(p, e, z) and lF(p, e, z) are governed by 

iv F p z d z T = -2eT + A(p, z) + JA(p, z) f J", (2) 
ivppAF = 2ef - A{p, «)t - JA(p, z) J", (3) 

where Up is the Fermi velocity and A(p, z) is a gap ma- 
trix. The quasiclassical Green's functions have symme- 
tries 



g(p,e,z) = g(p,e*,z)\ 
f(p,e,z) = -f(-p,-e,z) T . 



(4) 
(5) 



Here the superscript T denotes matrix transpose. Equa- 
tions Q and ([5]) give the general relation between the 
retarded Green's function (e = E + iO) and the advanced 
one (e = E—iO), and also between the Matsubara Green's 
functions on the positive and negative imaginary axes in 
the complex e plane. 

From g, one can calculate the local density of states, 



n(p, E, z) = Im 



-Trg(p,E + iO,z) 



(6) 



The even- and odd-frequency pair amplitudes are defined 
in terms of / by 



f^(p,e,z) 11 



r Y {p,e,z) 



-e,*)]. (7) 



In the following analysis, we consider the unitary states 
in which AA^ = |A| 2 is proportional to the unit ma- 
trix. The unitary condition is satisfied for all of the spin- 
singlet states and for many of spin-triplet states, e.g., 
those realized in the A and B phases of superfluid 3 He 4 
and the two-dimensional chiral p-wave state discussed for 
Sr 2 Ru0 4 £ 

In the unitary states, T and T for z — > oo take the 
form 



A°°(p) 
+ ifi(p,e) 



A°°(p)t 



(8) 



Here A°°(p) = A(p,z — > oo) is the bulk gap matrix 
and fl(p,e) = (|A°°(p)| 2 - e 2 ) 1/2 in which \A°°(p)\ corre- 
sponds to the bulk energy gap. The J 700 and J 700 make 



the right-hand side of Eqs. ([2]) and (J3)) in the bulk re- 
gion vanish and substitution of them into Eq. (fTJ) gives 
the well-known bulk solution of the quasiclassical Green's 
functions, g°° = e/fi(p,e) and f°° = A°°(p)/Cl(p,e). 

The boundary problem in the semi-infinite system can 
be solved in the following way^ To obtain the incom- 
ing (outgoing) solution J-p r< o (-7 r p^>o) J the differential 
equation @ is integrated from right (left) to left (right). 
Equation ([3]) for T is integrated in the opposite direction. 
The information on the surface is therefore included in 
the initial values of the outgoing solutions (J r p z >o and 
J r p i <o), while the incoming ones (J r p i <o and -Fp z >o) & re 
independent of the boundary condition at the surface 
apart from the implicit dependence through A(p, z). 



III. RELATION BETWEEN THE PAIR 
AMPLITUDE AND THE LOCAL DENSITY OF 
STATES 

We now consider the retarded Green's functions (e = 
E + iO) for the spin-triplet unitary states and show that 
the odd-frequency pair amplitude is closely related to the 
local density of states in the low energy limit. The spin- 
triplet states are characterized by the gap matrix of the 
form 



A(p,z) = d(p,z) ■ eria 2 , 



(9) 



where <x = (cr±, 02, 03) is the Pauli matrix. The d- vector 
in the unitary states can be expressed as the product of 
a real vector and a phase factor: 



d(p,z) = d ( P ,zy^- 



(10) 



In the absence of the magnetic field, the matrix J- does 
not have a spin-singlet component and we may put 



T = T ■ cri<T2, J~ - 



(11) 



It is straightforward to show that the incoming solutions 
of T and T for \E\ < |A°°(p)| can be parametrized as 



"33 z <0 



p 2 >o = «ae 



(12) 



where a is a unit vector and 7 is a real function; those 
are determined from 

v F p z d z d = 2 cos(7 - ipo) [do - (do ■ a)d] , (13) 
VFP z d z j = 2 [E - sin(7 - <^o)(cio • a)] , (14) 



and the initial values in the bulk region, 

a oo = -sgn(p 2 )d oo /|d oo |, 

7 00 = ~ sgn( Pz ) arcsin (E/\d™ 



(15) 
(16) 



Equations ([TJ and (JTTJ) yield g — go+g-cr and / = f-cri<j2 
with 



(17) 
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and 



T - T T 



(18) 



/e ZT ) . (20) 

(p,E+iO,z) 



where N = 1 — 2^" ■ T + T T . From the above expres- 
sions for g and /, one can derive the following relations: 

.90^ -i(/-^ 2 /), .90^ =-i(/-^ 2 /)- (19) 
Combining Eqs. (fT2j) and (fT9|) . we obtain 

(3oa)(p,B+io, z ) = 2 (/ e ~' 
Taking the imaginary part of this equation, we get 

n(p,E,z)\ {E{<Aoomi = \Re[D(p,E,z)e-^]\ (21) 
with 



D(p, E,z) = - f(p, E + iO, z) + f(p, E + iO, z) 



= - [f(p, E + iO, z) - f(p, E - iO, z)] , (22) 

where the last equality follows from Eq. ([5]). At E = 0, 
the right-hand side of Eq. (|2Tj) can be written in terms 
of the odd-frequency pair amplitude, 



n(p,0,z)= |Re[/ OF (p,iO,z)e-^] 



(23) 



n(p,0,z) = \Ref Ub {p,iO,z) 



It follows that the observation of a finite density of states 
at zero energy gives a direct evidence of the existence of 
the odd-frequency Cooper pairs. 

Note that the phase 7 at E = arises essentially due to 
the breaking of time-reversal symmetry. This is because 
the d- vector for the time-reversal invariant states can be 
taken to be real (tpo = 0) and then we see from Eqs. (p~4)) 
and (fl~6|) that j(E = 0) = at any position z. Thus, for 
the time-reversal invariant states, Eq. (|2"31 is reduced to 



(24) 



Equation (f24| can be applied to superfiuid 3 He-B, as dis- 
cussed below. 

Before proceeding to discussion on superfiuid 3 He, we 
briefly mention the relations similar to the above in the 
case of spin-singlet states with A(p, z) = Ao(p, z)ioi [see 
also Ref. in which an analog of Eq. ([20)) for a spin- 
singlet superconductor can be found] . For \E\ < |Ag°(p)|, 
we can put .Fp 2 <o = —ise^io-i and J r p z >o = ise~ 1 ~ i io2 
with s — sgn[Re A (p, z)]. Then, we can arrive, in a 
similar way to the spin-triplet 



n(p,0,z) 



|Re[/ OF (p,z0,z)e-^] 



(25) 



where /(P F is the odd-frequency spin-singlet pair ampli- 
tude, i.e., / OF = /o >F itT 2 . For the time-reversal invariant 
states, we can again take j(E = 0) to be zero. 
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FIG. 1. (Color online) Spatial dependence of |/j_,;(_E + iO, z)\ 
at E = 0.5Ab in the semi-infinite superfiuid 3 He-B at temper- 
ature T = 0.2T C (T c : superfiuid transition temperature). The 
surface specularity S is taken to be (a) S = 1 (specular limit) 
and (b) S — (diffusive limit). Self-consistently determined 
gap functions A\\ t ±(z) are used in the calculations. 



Let us turn to the surface odd-frequency state in super- 
fluid 3 He-B. The <i-vector in the semi-infinite superfiuid 
3 He-B has the forr 1 ^ 



d= (An(z)p x ,A l i(z)py,A±(z)p z ). 



(26) 



Here Ay t ±(z) are spatially dependent gap functions, both 
of which tend for z — > 00 to the isotropic bulk gap of the 
B phase, A°°(p) = A^. Since Aim (z) can be taken 
to be real, this state is an example of the time-reversal 
invariant spin-triplet states. 

The p-wave gap functions Aim (z) have the following 
spatial structure depending on the boundary condition 
at the surface?""— When the surface is specular, the 
perpendicular component A±(z) is strongly suppressed 
near the surface by the destructive interference effect due 
to quasiparticle reflection, while the parallel component 
A || (z) is slightly enhanced to compensate the condensa- 
tion energy lost by the suppression of Aj_(z). At rough 
surface causing diffuse quasiparticle scattering, the paral- 
lel component also suffers from the interference effect, so 
that both of Ay (2) and Aj_(z) are suppressed. The nu- 
merical results for An ±(z) determined self-consistently 
from the gap equation can be found in Refs. [T3 - [l^ . 

The self-consistent solution of / takes the form 



/ = (/|| cos 4>, /|| sin 4>, f± 



(27) 



where 4> = &rct&a{p y / p x ) is the azimuthal angle of p. 
The pair amplitudes /mi are a function of (p z ,e,z) and 
are independent of <p because of the rotational symmetry 
around the surface normal. 

To illustrate the existence of the surface odd-frequency 
Cooper pairs in superfiuid 3 He-B, let us consider the pair 
amplitude f±(p z , e, z) and its partial- wave components 



f±,l( e > z ) = {f±(Pz,£,z)Pi(p z ))p-, 



(28) 



where Pi(p z ) is the Legendre polynomial and {■■■)§ de- 
notes the angle average over the Fermi surface. It fol- 
lows from the general symmetry relation §5§ that the 
even-i (odd-i) components have the symmetry of odd 
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E/A B 

FIG. 2. (Color online) |Re/ OF (p, E + iO, z = 0)| in the semi- 
infinite superfluid 3 He-B at T = 0.2T C as a function of E/Ab- 
The inset is the surface density of states, n(p, E, z — 0). The 
two quantities are independent of the azimuthal angle (j> — 
arctan(p y / p x ) because of the rotational symmetry around the 
surface normal. The polar angle 9 — arccos(p z ) is taken to 
be 7r/6. The numerical results are shown for several values of 
the surface specularity S. 
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FIG. 3. \Ref OF (p,E + i0,z = 0)| shown in Fig. [2] for the 
diffusive limit (S = 0) is compared with the surface density 
of states, n(p, E, z — 0). 



(even) frequency. In Fig. [TJ we plot the magnitude of 
the self-consistent solution of f±,i{E + i0,z) at E = 
0.5 As as a function of z scaled by the coherence length 
Co = vf/2ttT c . The rough surface effect is taken into ac- 
count using the random 5-matrix theor y 17 ' 27 The sur- 
face roughness is parametrized by specularity S (0 < 
S < 1; S = 1 corresponds to the specular surface and 
S — » to the diffusive limit where the quasiparticlcs 
are scattered isotropically) i 20 i 21 The odd- frequency even- 
l Cooper pairs have substantial amplitudes at the surface 
in both of the specular (S = 1) and diffusive (S = 0) lim- 



its. At specular surface, in particular, f± has only the 
odd- frequency components. This is due to the reflection 
symmetry fx(ftz> e, 0) = /j_(~p 2 , e, 0) at the specular sur- 
face. 

In Figs. [2] and [3 the E dependence of |Re/ OF | for e = 
E+iO at the surface of superfluid 3 He-B is shown together 
with the surface density of states (SDOS), n(p, E, 0). The 
two quantities have quite similar E dependence. In the 
limit E — > 0, the both coincide exactly with each other, 
as we have shown analytically [see Eq. (|24p]. 

In Fig. the sharp peak (solid line) below the bulk 
gap (E/Ab < 1) corresponds to the Andreev bound 
state formed near the specular surface (S = 1). As 
the specularity S decreases, the peak is broadened and 
the density of states around zero energy increases. Such 
low energy excitations in superfluid 3 He-B have re- 
cently been observed by transverse acoustic impedance 
measurementsj 20 ' 21 In the experiment, the specularity 
is controlled by coating the surface with thin 4 He lay- 
ers. The specularity dependence of the surface midgap 
density of states is detected as a distinctive change in a 
peak structure found in the frequency dependence of the 
impedance below the bulk gap Ag. 18 ' 20 ' 21 This experi- 
ment gives a strong evidence of not only the formation of 
the surface Andreev bound states but also the existence 
of the surface odd-frequency Cooper pairs. 



IV. CONCLUSION 

Based on the experimental result a 20 ' 21 of the transverse 
acoustic impedance in superfluid 3 He and on the formula 
([24]) relating the odd-frequency pair amplitude to the 
zero-energy density of states, we have demonstrated that 
the odd-frequency Cooper pairs are formed near the sur- 
face of superfluid 3 He-B. The surface odd-frequency state 
is the origin of the anomalous proximity effect predicted 
for a diffusive normal metal/spin-triplet superconductor 
junction^ and an analogous system composed of aerogel 
and superfluid 3 He-B^ Since the pairing symmetry of 
bulk superfluid 3 He is well established, its surface or in- 
terface provides an ideal environment for studying the 
physics of the odd-frequency state. 
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